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OA ■ Abstract. Bose-Einstein correlations of identical charged-pion pairs produced in hadronic Z decays 

£H ' are analyzed in terms of various parametrizations. A good description is achieved using Levy stable 

distributions. The source function is reconstructed with the help of the T-model. 
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INTRODUCTION 



In particle and nuclear physics intensity interferometry provides a direct experimental 
method for the determination of sizes, shapes and lifetimes of particle-emitting sources 
(for recent reviews see JH0]). In particular, boson interferometry provides a powerful 
tool for the investigation of the space-time structure of particle production processes, 
since Bose-Einstein correlations (BEC) of two identical bosons reflect both geometrical 
• and dynamical properties of the particle radiating source. 

In e + e annihilation BEC are maximal if the invariant momentum difference is small, 
even when one of the relative momentum components is large, as was seen by TASSO 
J3I] and which we have confirmed. For a hydrodynamical type of source, on the contrary, 
BEC decrease when any of the relative momentum components is large JUS]. 

Here we investigate various parametrizations and find that a good description of 
the Bose-Einstein correlation function can be achieved using Levy stable distributions 
as the source function. Within the framework of models assuming strongly correlated 
coordinate and momentum space, we then reconstruct the complete space-time picture 
of the particle emitting source in hadronic Z decay. 

For our analysis we use a sample of about 500 thousand two-jet events, selected by 
the Durham algorithm o] with y cut = 0.006, from e + e annihilation data collected by 
L3 at a center-of-mass energy of 91.2 GeV. 

PARAMETRIZATIONS OF BEC 

The two-particle correlation function is defined as: 



(1) 



where pi(p\,P2) is the two-particle invariant momentum distribution, pi (pi) the single- 
particle invariant momentum distributions and pi the four- momentum of particle i. Since 
we are only interested in BEC, the product of single particles densities is replaced by the 
so-called reference sample, po(pi,P2), the two-particle density that would occur in the 
absence of Bose-Einstein interference. Here we use mixed events as a reference sample. 

After some assumptions yj,|2|], this two-particle correlation function is related to the 
Fourier transformed source distribution. In this case 

R2(puP2) = 1 + \HQ)\ 2 , (2) 

where f(x) is the density distribution of the source, Q is the invariant four-momentum 
difference, Q = —(pi — pi) 2 and f(Q) is the Fourier transform of f(x). 



Gaussian distributed source 

The simplest assumption is that the source has a symmetric Gaussian distribution, in 
which case f(Q) = exp {ip,Q — J an ^ 

Ri(Q) = y [l + A exp {-(RQ) 2 )} (1 + SQ), (3) 

where the parameter 7 is a constant of normalization, A an incoherence factor, which 
measures the strength of the correlation, and (1 + 5(2) is introduced to parametrize 
possible long-range correlations not adequately accounted for in the reference sample. 

A fit of Eq.© to the data results in an unacceptably low confidence level from which 
we conclude that the shape of the source deviates from a Gaussian. The fit is particularly 
bad at low Q values. 



Levy distributed source 

Adopting Nolan's S(a,/3 = 0,y, 8; 1) convention [6] for the symmetric Levy stable 
distribution with rescaling of the scale parameter y to R and the location parameter 5 to 
xo, the Fourier transform (characteristic function) f(Q) has the following general form: 

f(Q)=exp(iQx -\RQ\ a ). (4) 

The index of stability, a, satisfies the inequality < a < 2. The case a = 2 corresponds 
to a Gaussian source distribution. For more details seeja]. 
Then R2 has the following, relatively simple form [7]: 

^2(2) = y[l + Aexp(-( J R<2) a )](l + 5<2). (5) 

After fitting Eq.© to the data it is clear that the correlation function is far from Gaussian: 
a ~ 1.3. The confidence level, although improved compared to the fit of Eq.©, is still 
unacceptably low. 
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FIGURE 1. The Bose-Einstein correlation function R2- The curve corresponds to the fit of Eq.(|6jl. 

Since there is no particle production before the onset of the collision, a more appro- 
priate form of the source distribution for the time component is the asymmetric stable 
distribution. In this case, one obtains the following result for the correlation function 0]: 

Ri(Q) = 7 [1 + A cos l(R a Q) a ] exp (- (RQ) a )] (1 + 8Q), (6) 

where R a is an additional parameter, a measure of the onset of particle production 

The fit of Eq.® to the data, shown in Figure[H is statistically acceptable. The data are 
well described by the fit. Note that for Q between 0.5 GeV and 1 .5 GeV the data points 
go below the dashed line, which stands for the long-range correlations extrapolated to 
lower Q values. These data points indicate an anti-correlation in the Q 1 GeV region. 
This property of the data is well reproduced by the fitted curve, which also goes below 
unity as a result of the cosine term in Eq.©, which comes from the asymmetric Levy 
assumption. The fitted value ofais0.82±0.03. 




THE t-MODEL 

A model of strongly correlated phase-space was developed [9] to explain the experimen- 
tally found invariant relative momentum Q dependence of Bose-Einstein correlations in 
e + e~ reactions. This model also predicts a specific transverse mass dependence of R2, 
that we subject to an experimental test here. 

In this model, it is assumed that the average production point x^ of particles with a 
given momentum k* 1 is given by 

x^(k^)=dk^. (7) 
In the case of two-jet events, d = where % = \/t 2 — k\ is the longitudinal proper- 
time and mj = \ m 1 + p\ is the transverse mass. The second assumption is that the 



distribution of x^(k^) about its average, 5a (x^(k^) — x^(k^)), is narrower than the 
proper-time distribution. Then the emission function of the T-model is 

p CO 

S(x,k) = / dTH(T)8b(x-dk)Ni(k), (8) 
Jo 

where H(t) is the longitudinal proper-time distribution, the factor 5a (x — dk) describes 
the strength of the correlations between coordinate space and momentum space variables 
and A^i (k) is the experimentaly measurable single-particle spectrum. 

In the plane-wave approximation the Yano-Koonin formula fioll gives the following 
two-pion multiplicity distribution: 

Pl{k\M) = J d 4 Xid 4 X2S(xi,ki)S(x2,k 2 ) (1 +cos[(fci -k 2 )(xi -x 2 )]) . (9) 

Approximating the 5a function by a Dirac delta function, the argument of the cosine 
becomes 

(*l -k 2 )(x l -x 2 ) = -0.5 (di +d 2 )Q 2 . (10) 
Then the two-particle Bose-Einstein correlation function is approximated by 

RiihM) = 1 +ARe£ 2 (JL\ , (11) 

where H(co) = J dxH(z) exp(z'coT) is the Fourier transform of H[x). Thus an invariant 
relative momentum dependent BEC appears. 

Guided by the result of the previous section, we use an asymmetric Levy distribution 
for the longitudinal proper-time density. Thus the corresponding BEC function has an 
analytic, although somewhat complicted form ||7l|8|]: 



R2(Q 2 ,m T ) =7 



1 + A cos — YA — — — exp 



B (12) 



where the parameter To is the proper-time of the onset of particle production, At is a 
measure of the width of the proper-time distribution, A = tan f^pj and B = (1 + 8Q) 

After fitting for various mj interval we find that the quality of the fits is statistically 
acceptable and the fitted values of the model parameters are stable and within errors 
the same in all investigated mj interval. The T-model with a one-sided Levy proper- 
time distribution describes the data with parameters To = fm, a ss 0.80 ±0.05 and 
AT^2.0±0.6fm. 



RECONSTRUCTION OF THE EMISSION FUNCTION 

In order to reconstruct the space-time picture of the emitting process we assume that the 
emission function can be factorized in the following way: 

S(r,z,t)=I(r)G(T])H(T), (13) 



FIGURE 2. Two views of the longitudinal part of the source function normalized to the average number 
of pions per event. 



where /(r) is the single-particle transverse distribution, G(r\) is the space-time rapid- 
ity distribution of particle production, which approximately coincides with the single- 
particle rapidity distribution, and H(t) is the observed proper-time distribution. 

With these assumptions one can reconstruct the longitudinal part of the emission 
function integrated over the transverse distribution. It is plotted as a function of t and 
z in Figure El It exhibits the typical boomerang shape with a maximum at low t and z 
but with tails reaching out to very large t and z values, a feature already observed for 
hadron-hadron 111 ill and heavy ion collisions lll2ll . 

The transverse profile, which follows from Eq. ©, has the following form: 



This equation describes the particle production in coordinate space as a function of the 
proper-time T. It describes the expansion of the source as the proper-time increases. 
The particle production probability is proportional to the proper- time distribution H{x). 
Figure |3] shows the transverse part of the emission function for various proper-times. 
Particle production starts immediately, increases rapidly and decreases slowly. A ring- 
like structure, similar to the expanding, ring-like wave created by a pebble in a pond, 
is reconstructed from L3 data, as shown in Fig. 3. An animated gif file that shows this 
effect is available from H 1 311 . 




(14) 
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FIGURE 3. The source function normalized to the average number of pions per event for various 
proper- times. 
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